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A WALL-WAKE VELOCITY PROFILE FOR TURBULENT 
COMPRESSIBLE BOUNDARY LAYERS WITH HEAT 
TfWNSFER 

I. INTRODUCTION 

The wall -wake velocity profile for isoenergetic compressible turbu- 
lent boundary layer flow developed by Mathews et al., [1], and modified by 
the present authors, [2], was based on the work of Maise and McDonald [3]. 
In their work, Maise and McDonald suggested that the defect of Van Driest's 
[4] generalized velocity u* has the same form for compressible turbulent 
boundary layer flow with zero pressure gradient as the velocity defect in 
incompressible flow given by Coles [51. Mathews et al., [1], went further 
to suggest that a boundary layer profile similar in form to Coles' wall- 
wake profile, but the velocity replaced by u*, might be applicable for com- 
pressible flow with pressure gradient. The resulting profile and the sub- 
sequently modified form of it have indeed been round to provide a good re- 
presentation of the mean velocity distribution for adiabatic flows with 
pressure gradient. 

Since Van Driest's generalized velocity u* can take into account the 
effect of heat transfer, it was reasonable to investigate the use of pro- 
files expressed in terms of u* for non-adiabatic flows, f'aise and !*cDonald 
[31 attempted a u*-defect correlation for such flows using the seme rela- 
tionship they had used for adiabatic flat-plate flows. They examined a 
number of non-adiabatic velocity profiles but found an unacceptable level 
of scatter in the data. 

Lewis, Gran and Kubota [6], working with Van Driest's transformed 
velocity, U*, which differs from u* by a constant for a given profile. 
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were able to correlate low-speed and compressible turbulent boundary layer 
data for adiabatic flow with p '-assure gradients through the use of a pres- 
sure gradient parameter - (fi*|^/T^)(dp/dx) suggested by Alber and Coats 
[7]. Subsequently Gran et al . , [8] used the wall-wake profile proposed in 
Reference 6 to represent boundary layer profiles for cold-wall flows sub- 
jected to pressure gradients. They found that the pressure gradient para- 
meter correlated low-speed data with both adiabatic and cold-wall high- 
speed data. 

The velocity profile used in References 6 and 8 has a non-zero velo- 
city gradient at the boundary layer edge. The purpose of this note is to 
suggest a modified form of the wall -wake profile which is applicable to 
flows with heat transfer, and for which 3u/3y = 0 at y = o. The modified 
profile also takes into account the effect of turbulent Prandtl number. It 
has be*- found to provide a good representation of the experimental data 


fror everal sources and the values which are determined by a least 
sqjdres fit of the profile to the data agree well with values which were 
measured by the floating element technique. 


II. VELCCITY PROFILE 


For incompressible flow the wall-v/ake velocity profile as prcoosed 
by Coles [5] may be written as, 

u/u^ = (1/K) In (yu^/v) + C + nW(y/o)/K (1) 

Setting u = u^ and V.'(y/6) = 2 at y/<5 = 1 in Eq. (1) and subtr'acting 
the resulting equation from Eq. (1) leads to an expression for the velocity 
defect, 

(Ug - u)/u^ = - (1/K) In (y/6) + (n/K) [E - l'(y/0] 


( 2 ) 
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Maise and McDonald [3] suggested that for compressible flow the 
generalized velocity defect also assumes a unique universal functional 
relationship such as that given by Eq. (2), i.e., 

(u* - u*)/u^ = - (1/K) In (y/6) + (n/K) [2 - W(y/«;)] (3) 

They found that this expression provided a gjod repr»“sentation of adia- 
batic, flat-plat'i, boundary layer profiles with n/K set equal to 1.25. 

Mathews vX al., suggested for isotnergetic flow a profile similar 
in form to that given by Eq. (1) but witi u replaced by u*, i.e., 

u*/u^ = (1/K) In (yu^/v^) + C + (ii/K) W(y/6) (4) 


This profile was found to be applicable to adiabatic flows with pressure 
gradient for a range of Mach numbers and flow configurations. The value 
of n/K, which was determined through a least squares fit of the wall-wake 
profile to experimental d'ta, was different for different values of pres- 


sure gradient. 
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layer edge has a non-zero value. This shortcoming is avoided with the 
modified profile proposed by Sun and Childs [2]: 
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With the value of a taken to be 1. 

Equation (5) has been tested against a large number of experimental 
boundary layer profiles on adiabatic walls for a wide range of Mach 
numbers and has been found to provide a good representation of the velo- 
city profile. However, for flow on a non-adiabatic wall, the applicability 
of Eq. (5) deteriorates as the ratio of the temperature on the wall to the 
free stream total temperature decreases. 

Van Driest [5] developed a compressible law-of-the-wall from the 
concept of mixing length in a manner similar to that followed by Prandtl 
for incompressible flow. It may be written as, 
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A comparison of Eq. (8) with Eq. (1) suggests that a boundary layer profile 
for compressible flow with heat transfer should be written as. 
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The left-hand side of this expression is the transfcrmied velocity, U*, 
which is used in References 6 and 8. 

In Van Driest' s development of Eq. (8), he assumed a turbulent Prandtl 
number of unity, which leads to the following temperature and density dis- 
tributions across the boundary layer, 
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In addition, he assumed a constant shear stress near the wall. Thus, 
with Eq. (9) as written, the velocity gradient at the boundary layer edge 
will have a non-zero value. If we assume, as In Reference 2, that t may be 
expressed as t = t (1-n^), and If for Pr^ M » the temperature and density 
distribution': may be written as (cf. Schllchting [9] ), 
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we may follow the same procedure as is used in Reference 2 to obtain a 
boundary layer velocity profile of the following form, 
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For a which corresponds to the assumption of a linear shear stress 
distribution in the derivation of the compressible law-of-the-wall , and 
for Pr^ 1, u** u* as in Eq. (9). Setting n * 1 in Eq. (14) yields 
the following expression for n/K, 
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For mathematical convenience [1-3] we may replace [2-W(n)l in 
Eq. (15) by [1 + cos nir] and write the wall -wake velocity profile as, 

U 2A2-B, 

iT * — JT'z s1n - l1 

X !U ^ (In r, + I - I In (l+O-r,®)’*)) 


n “t h 
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(17) 


with n/K given by Eq. (16). A least squares fit of the profile to experi 
mental data yields values of 6 and u^ or C^. As Pr^ 1, u** u*, so 

that Eq. (17) reduces to Eq. (5). 

The value of n/K may be calculated from the conditions at the v/all 
and at the boundary layer edge by Eq, (16). However, it may also be 
obtained from velocity profile data by integrating Eq. (15) from n = 0 to 
n = 1. For a = 1 the integration yields, 
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III, DISCUSSION OF RESULTS 

The method of let'st squares has been used to fit the wall -wake pro- 
file, Eq. (17), to experimental velocity profiles reported by Hopkins and 
Keener [10], Voisinet et al., [11, 12], Horstman and Owen [13] and Kilburg 
[14] for zero pressu»*e gradient flows with heat transfer. The wall -wake 
velocities were obtained under the assumption of Pr^ * 0.8 and a « 1. 

K was taken as 0.4 and C as 5.1. Data in the sublayer were excluded. The 
results shown here are from studies for which directly measured values of 
wall shear stress were reported [10, 11, 13]. The data reported by 
Voisinet et al., and Hopkins and Keener were obtained on wind tunnel walls. 
The Horstman-Owen data are for flow along the cylindrical surface of a 
cone-cylinder model whose axis was aligned with the primary wind tunnel 
flow direction. 

Comparisons of the experimental and wall -wake profiles are given in 
Figure 1. Also shovm on the figures are values of 6*, e, and n/K. 

The wall -wake profile is seen to provide a good representation of the 
experimental velocity profiles. In addition, +he values determined by 
the fit of the proposed profile agree reasonably well with the directly 
measured values. For reasons which are not known at this time the wall- 
wake profile does not provide as good an overall fit to the data reported 
by Voisinet et al., as it does to the data of Hopkins and Keener or Horst- 
man and Owen. As Figure 1 shows, the wall-wake profile agrees quite well 
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V 

with the Volsinet data bj^twcen the wall and the value of 6 determined 

* 

by the least squares fit. However, at the value of 6 determined by the 

fit, u/Ug * 0.965 rather than 0.995 which Is sometimes t^ken to correspond 

to the boundary layer edge In high speed flows. For the «.4,92 flow 

examined by Volsinet, P^/Pq ■ 0.56 and T^/T^ ■ 0.938 at u/u^ ■ (J;96J. 

e e 

On the other hand, for the Hopkins-Keener and the Horstman*0wen datd the 

. f 

values of u/u^, Pq/Pq * ^o^^o essentially unity at the value of 6 

e e 

determined by the fit. Values of 6 were also determined for a • In 

Eq. (17). As would be expected, they are lower than for a ■ 1 and u/Ug, 

^o^^o ^o^^o ^ correspondingly lower. 

*^e e 

The values of n/K shown In Figure 1 were determined from Eq. (16), 
assuming Pr^ * 0.8 and a * 1 and with u^ and 6 determined from the wall- 
wake fit. Although the pressure gradient was zero for the flows under con- 
sideration, the value of n/K varies considerably from profile to profile 
and Is generally higher than for low-speed or compressible adiabatic flows. 
Note that, for the modified profile used here, the expression for n/K 
differs from that used in References 1, 3, 6 and 8 through the Inclusion 
of the Prandtl number effect, through the term 0.614/aK, and through the 
differences in and 6 which are obtained with the modified profile. 

Figure 2 shows a comparison of the experimental data from References 
10, 11 and 13 with the u**-defect relationship given by Eq. (15) for 
Pr^ * 0.8 and with the value of K taken as 0.4. The values of r./K which 
ere listed were calculated from Eq. (18) using the measured values of 
shear stress to evaluate u^. The experimental values of u** have been 
calculated from the reported values of Mach number and total temperature. 

As is shown, Eq. (15) provides a reasonably good representation of the 
u**-defect for each profile. 
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The results may also be plotted In the form (u** - u**)/u vs. 
y/A, where a Is analogous to the defect thickness used by Gran et al., [8] 

In their comparison of cold-wall data with low-speed data [15]. Such a 
plot is shown in Figure 3. It is seen thU the data of Hopkins and Keener 
agree quite closely with the low-speed results. The Horstman-Owen data* 
for which n/K is about thirty percent higher than for the Hopkins Keener 
profile, deviates from the low-speed results in the outer part of the 
boundary layer. The Voisinet data, which n/K is considerably higher than 
for either of the other two profiles, are quite different from the low- 
speed data, even though the wall -wake profile fits the experimentally de- 
termined velocity distribution quite closely between the wall and the value 
of 6 determined by the least squares fit. 

The values of 6* and e which were determined from the wall-wake pro- 
file are listed on Figure 1, along with the experimental values reported 
in the references. The results agree quite well for the Horstman-Owen pro- 
file, but only moderately well for the other two. In the i-ase of the data 
by Voisinet, the difference appears to result from the fact that the velo- 
city, total temperature and, total pressure continue to increase well beyond 
the value of 6 determined by the least squares fit. For the data of 
Hopkins and Keener on the other hand, the difference is due to differences 
between the experimental density profile and the profile implied by Eq. 

(13). Data reported by Schlichting [16] suggests that the turbulent Prandtl 
number for high-speed flow on a flat plate is on the order of 0.9 and that 
for flow on the wall of a plane channel is on the order of 0.8. Recent 
work of Heier et al., [17] also suggests that Pr^ is near 0.8 for flow on 
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thf wall of a plane channel. The total temperature o.strlbutlon reported 
by Hopkins and Keener Imply a considerably lower Prandtl number, and con- 
sequently, a considerably different density distribution than Is obtained 
when Pr^ ■ 0.8 Is assumed In Eq. (13). This In turn leads to differences 
between the experimental and the wall -wake values of 6* and e. 

IV. CONCLUSION 

An extension of the wall -wake velocity profile of Reference 2 to 
compressible turbulent boundary layer flow on non-ad1abat1c walls has been 
tested against experimental data covering a wide range of Mach number and 
heat trarsfrr and has bee-^ found to provide a good representation of experi- 
mental velocity profiles. The profile considers the effect or turbulent 
Prandtl number and allows for variation of the shear stress near the wall. 
Values of C^ determined by a least squares fit of the .vall-wake profile to 
experimental data agree reasonably well with directly measured values. 

The values of ^ determined by the fit are closer to those based on u/u^ » 
0.995 than are found with earlier versions of the v/all-wake profile which 
were developed under the assumption of a constant shear stress near the 
wall. The values of the integral properties 6* and g as determined with 
the wall-v/ake profile deviate miore, In some Instances, from experimentally 
determined values than was found to be the case for adiabatic flows. This 
occurs because the density profiles given by Eq. (13) and those computed 
from the experimental total temperature and Kach number profiles differ 
considerably. 


n 


The wake function coefficient n/K, varies considerably more than 
has been reported for zero-pressure-gradient non-adiabatic flows in 
earlier references, and* for most of the data examined, is higher than 
has been reported for adiabatic flows. 
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Fig. 1 

Fig. 2 
Fig. 3 


Velocity distributions in boundary layer with heat 
transfer. 


Defect of generalized velocity profiles. 


Comparison of generalized velocity defect for low- 
speed flow and compressible flew with heat transfer. 
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